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ABSTRACT 

We explored a method to reconstruct the distribution function of the Galac¬ 
tic thick disc within the action space where nearby thick-disc stars are distributed. 
By applying this method to 127 chemically-selected thick-disc stars in the Solar 
neighbourhood, we found that the vertical velocity dispersion that corresponds to 
the reconstructed distribution function declines approximately as exp(—at 
4.5 kpc < i? < 9.5 kpc, with = 8.3 ± l.l(rand.) ± 1.6(sys.) kpc. Also, we found that 
the vertical velocity dispersion tjz of our local thick-disc stars shows only weak depen¬ 
dency on radial and azimuthal velocities [vr,v^). We discuss possible implications of 
these results on the global structure of the Milky Way thick disc. 


Key words: - Galaxy: 
kinematics and dynamics 


1 INTRODUCTION 

Historically, the thick disc of the Milky Way was first iden¬ 
tified through star counts toward the Galactic Poles as a 
verti cally extended disc component wit h scale height ~ 1 
kpc (lYoshiill 19821 : iGilmore fc Reidlll983h . in addition to the 
thin disc with scale height ~ 0.3 kpc that dominates the 
disc stars in the immediate Solar neighbourhood. Therefore, 
at the outset, the tentative definition of the thick-disc stars 
was those disc stars with large vertical orbital excursions 

and large vertical velocities. _ 

Later, spectroscopi c studies fe.g.. I^nsbv et all l2003l . 
l2004l : lReddv et 51l2006l i on kinematically-selected disc stars 
suggested that stars with large vertical motions (which are 
likely to belong to the thick disc) tend to show lower [Fe/H] 
and higher [a/Fe] than those with small vertical motions 
(thin-disc stars). These chemical properties suggest that the 
thick-disc stars are older than thin-disc stars (lower [Fe/H]) 
and that the star formation timescale of the thick disc was 

shorter than that of the thin disc ( higher [g/Fe]). _ 

Recently, Klaus Fuhrmann jFuhrmannl 1 19981 . |2004 
I 2 OO 8 I . [2oT3) investigated a kinematically-unbiased volume- 
complete sample of Solar-type disc and halo stars located 
within 25 pc of the Sun. The distribution of his sample 
stars in the [Fe/H]-[Mg/Fe] space reveals two che mically dis¬ 
tinct populations of disc stars (see Figure 15 of iFuhrmannI 
H). Based on the kinematical properties of these popu¬ 
lations, he identified the lower-[Fe/H] and higher-[Mg/Fe] 
population to be the thick disc, and the other population 
to be the thin disc. This clear separation between the thin 
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and thick discs is also confirmed in the nearby (heliocen- 
tr ic distance d < 1 00 pc ) kinematically-unbiased sample 
of lAdibekvan et al.l ll2012l i for which [Fe/H] and [a/Fe] are 
available through high-resolution spectroscopy. These recent 
observations suggest that the thick disc is better defined 
by che mical compositions of stars, especially by [Fe/H] and 
[a/Fe] llMasserori k. Gilmordl2015h . 

In the past decades, our understanding of the struc¬ 
ture of the thick d isc towards the Galactic Poles has been 
greatly imp roved jGilmore et al.l 1 19891 : iRix fc Bovvl l2013l : 
lYoshii|[2013l L The next step forward is to unravel its more 
global structure, such as the radial dependence of its vertical 
structure. So far, many attempts have been made to fulfil 
this goal, and they are broadly categorised into two classes 
of studies. 

The first class of studies are based on nea rby samples 
of disc stars. One such example is iBinne 3 (120121 '). who fitted 
the distribution of local disc stars with his distribution func¬ 
tion model. He discussed some global structure of the thick 
disc by looking into his best-fit models. Although this kind 
of studies can predict almost everything if the best-fit mod¬ 
els are reasonably correct, one critical problem with these 
studies is the validity of the functional forms of the assumed 
distribution functions. 

The second class of studies are based on in-situ sam¬ 
ples of (relatively) distant disc stars. This class can be 
further categorised into three sub-classes: those studies us¬ 
ing (i) high-resolution spectroscopic samples; (ii) medium- 
resolution spectroscopic samples; or (iii) photometric sam¬ 
ples. The advantage of sub-class (i) studies is that we can 
define the thick disc purely by means of the stellar chemistry 
dBensbv et al.llioiH 'l. However, the number of stars that are 
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currently available is less than a few hundred, and this small 
sample size makes it hard to obtain some statistical proper¬ 
ties of distant disc stars. Also, since the errors in distance 
and proper motion are large, kinematical analyses are diffi¬ 
cult for these stars. In the sub-class (ii) studies, much larger 
sample of stars are available than in t he sub-class (i) stud - 
ies. A recent example of this sub-class is I Bow et al.l (l2012al l. 
who studied the density distribution of chemically defined 
disc populations by fitting the SDSS/SEGUE data with ana¬ 
lytic models of the density profile. However, since the chemi¬ 
cal abundances of these samples are less accurate than high- 
resolution samples, some blending of the thin- and thick-disc 
stars is inevitable in this kind of studies. Most of the sub- 
class (iii) studies are based on the star count method (e.g., 
iRobin et al.ll20l3 'l. Photometric samples are advantageous 
in grasping the global structure of the stellar disc because 
the available sample size is the largest among these sub¬ 
classes, and because the sample stars cover a wide range of 
heliocentric distance. However, since the photometric data 
lack chemical information for each sample star, it is not pos¬ 
sible to separate the sample into stellar components. This 
inseparability means that one needs to adequately model all 
the stellar components that contribute to the sample, which 
is harder than adequately modelling the thick disc only. 

By taking into account these problems, in this paper we 
explore the possibility of constraining the global structure of 
the thick disc based on a small number of nearby chemically- 
selected stars, but not relying on analytic modelling of the 
thick disc. Specifically, we reconstruct the distribution func¬ 
tion of the thick disc within a certain portion of the action 
space which is accessible from local observations. 

This paper is organised in the following manner. First, 
we describe our sample stars in section [2] In section |3] we 
introduce the concepts of the observable action space E0 and 
the observable distribution function /q. There, we discuss 
the possibility of inferring some information on the velocity 
dispersion of the thick disc by using /q. In section U] we 
present our method to reconstruct /q from a local sample 
of thick-disc stars. In section O we present our main results. 
The implications from our study are discussed in section [51 
and section 0 sums up our analyses. 


2 SAMPLE 

2.1 The coordinate system and some numerical 
constants 

In this paper, we assume that the Milky Way potential 
is axisymmetric and adopt the usual Galactocentric cylin¬ 
drical coordinate system {R, (j), z) to describe the three- 
dimensional positions of stars in the Milky Way. The three- 
dimensional velocities along the respective coordinates are 
denoted by {vfi,v^,Vz), where vr is positive outwards, 
is positive in the direction of Galactic rotation, and Vz is 
positive in the direction of the North Galactic Pole. 

We assume that the local standard of rest (LSR) is in 
a circular orbit w ith a rotation speed of ulsr = 240 km s“^ 
llReid et al.ll20lJ l and that the Galactocentric distance of 
the Sun is Rq = 8.3 kpc. We also adopt the peculiar mo- 
tion of the Sun with respect to the LSR determined by 
ISchonrich et al.l ll2010h . 


2.2 Construction of the local thick-disc sample 


In this paper, we assume that the thick disc of the Milky 
Way is a single stellar population and define the thick-disc 
stars to be those stars with —1.0 ^ [Fe/H] ^ —0.2 and 
[Mg/Fe] > 0.2, by taking into acc ount the Solar n eigh- 
bour observations of disc stars (e.g., lFuhrmannll201lf) . We 
constructed a kinematically-unbiased sample of 127 thick- 
disc stars, for which accurate position, velocity, [Fe/H], and 
[Mg/ Fe] are available. 84 st ars in our sample were taken 
from lAdibeky an et all (|2012|), a nd th e other 43 stars were 
taken from Soubiran fc Girai^ (1200Sfl . Most of our sample 
stars are distributed within 100 pc from the Sun. 


2.2.1 Sample from Adibekyan et al. (2012) 

lAdibekyan et all (I2OI2I') published a kinematically-unbiased 
sample of 1111 stars, for which accurate [Fe/H] and [Mg/H] 
estimated from high-resolution spectroscopic observations 
as well as the three dimensional Gartesian velocity based 
on Hipparcos observations are available. We derived the 
distances and the line-of-sight velocities for 1103 stars in 
their sample tha t were also included in the updated Hip¬ 
parcos catalog of Ivan LeeuwenI (l2007l l. In the course of this 
inversion, we took into account th e Solar peculiar motion 
adopted in lAdibekyan et al.l 1I2OI2II and w e used the three- 
dimen sional Cartesian velocity taken from I Adibekyan et al.l 
1I2OI2II and the position on t h e sky and the proper mo¬ 
tion taken from f^n Leeuw^ (l2007h . We then calculated 
the three-dimensional velocity {vr,v^,Vz) in our coordi¬ 
nate system, by using our assumption on (Ro,ulsr). From 
these 1103 stars, we selected those stars that satisfy both 
the thick-disc chemical criteria (—1.0 ^ [Fe/H] —0.2 
and [Mg/Fe] 0.2) and the kinematical criteria of Vcf, ^ 
70 km s“^ and ^Vr^Gj^ ^ 160 km s“^. The latter velocity 
cuts were applied to minimise possible contamination from 
halo stars. Finally, we discarded those stars that had been 
flagged as non-single stars, and we obtained 84 stars. 


2.2.2 Sample from Soubiran & Girard (2005) 

ISoubiran fc GirardI (1200 Sf) published a catalog of 743 stars 
compiled from 11 catalogs of nearby stars and listed the 
chemical abundances of these stars n Among these 11 
catalo gs, those cata lo gs tak en from lA llende Prieto et aP 


(2004ll. Chen et al.l (200^. lEdvardsson et al.l ( 1993lL 


iMishenina et al.l ( 2004 ), and Reddv et al.l ( 200; ) are 
kinematically unbiased, so we first selected those stars 
that are included in either of these five catalogs. Then, 
we excluded those stars that happened to be included 
in lAdibekyan et all (I2OI2I') . The position and velocity of 
these stars were then calculated eit her by using the data 
from Genova-Copenhagen Survey (iHolmberg et^aljjgOOgh 
or by using the updated Hipparcos catalog ( van LeeuwenI 
I2OO7II and SIMBAD database. Finally, as in the previous 
subsection, we applied the chemical and velocity criteria 
and obtained 43 thick-disc stars. 


^ When more than one published data are available for a certain 
chemical abundance of a given star, they list the averaged value. 
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Figure 1. Chemical abundances of our thick-disc sample. 


2.2.3 Caveats on the velocity cuts 

As described above, we adopted a pair of velocity cuts to 
exclude halo stars. Rigorously speaking, these cuts poten¬ 
tially could have introduced some kinematical bias on our 
analyses. However, we have confirmed that adopting a looser 
velocity constraint of > 0 km s“^ (without cuts on vr or 
Vz) has only negligible effects on our final results. 


2.3 Chemical and Kinematical properties of our 
sample 

The distribution of our sample stars in the [Fe/H]-[Mg/Fe] 
plane is shown in Figure [T] The uncertainties in the chem¬ 
ical abundances are small enough (typical uncertainties are 
A[Fe/H1 = 0.03 dex and A[Mg/Fe] = 0.05 dex for the 
lAdibekvan et al 1 l 2012 l sample) compared with the scatter 
in this diagram. Therefore, our sample does not suffer from 
severe contamination from low-[Mg/Fe] stars. 

Figure [ 2 ] shows the velocity distribution of our sample 
in the {v^,vr), (v^,Vz), {vr,Vz) and [v^, + vl) planes. 

Our sample stars show a mean velocity of {{vr), ivf), (Vz)) = 
(10.5±5.3, 192.3T3.6, —7.9T3.7) km s~^, and thus lag from 
the LSR by nearly 50 km s~^. The three-dimensional veloc¬ 
ity dispersion is (cTij, cr^) = (59.6,40.9,41.6) kms“^. 

In Figure [31 we show the value of Oz as functions of Vrj, 
and |uii|. To make these plots, we first sort the sample in 
V(f, or |uij|, evaluate the dispersion for a binned sample 
of A = 40 stars, and move through the sample in steps of 
4 stars (so that any pair of adjacent bins share 36 stars). 
The derived value of Oz for each bin is plotted against the 
median value of Vtf, or |ufl| of the bin. 

Although the bin size of 40 stars is rather large com¬ 
pared with the total sample size of 127 stars, the result of 
this figure suggests that there is no strong trend in az as 
a function of |u_r| or v^. This property suggests that the 
vertical kinematics (az) is not strongly correlated with the 
in-plane motion (|ufl| or v^). 
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Figure 3. Vertical velocity dispersion Oz of our thick-disc sample 
as a function of (a) and (b) |ufl|. Each data point represents 
the value of Oz of a binned sample of 40 stars, and each pair of 
adjacent bins share 36 stars in common. The uncertainty in Oz is 
calculated by assuming Gaussian distribution of vertical velocity 
Vz in each bin. 


3 THEORY 

3.1 Observable distribution function 


Most of the thick-disc stars in the Solar neighbourhood 
are on non-circular orbits. More specifically, the majority 
of nearby thick-disc stars have guiding radii smaller than 
the Galactocentri c distance of the Sun (R p), due to the 
asymmetric drift l|Binnev fc Tremainell200^ 'l. Therefore, if 
the thick disc is a static stellar system, kinematics of nearby 
thick-disc stars carry some information on the thick disc out¬ 
side the Solar neighbourhood, especially on the inner part 
of the thick disc. 

Of course, nearby thick-disc stars do not carry all the 
information of the thick disc, since these stars are inher¬ 
ently biased in that they represent a portion of the thick-disc 
stars whose orbits penetrate through the Solar neighbour¬ 
hood. In other words, a nearby sample only covers a cer¬ 
tain fraction of action space F© where the co rresponding or¬ 
bits b ring stars to the Solar neighbourhood (iMav fc BinnevI 
Il986h . Hence, the best we can do with nearby samples is to 
reconstruct the distribution function f(J) within this ob¬ 
servable action space r 0 , i.e. to reconstruct the observable 
distribution function defined by 


feiJ) 


f(J) (JgFo) 

0 (otherwise). 


( 1 ) 
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Figure 2. Velocity distributions of our thick-disc sample in (a) (u^,U 2 )-space, (b) (uj?,U 2 )-space, (c) (u^,'i;j 2 )-space, and (d) 
+^?)-space. 


We here note that /o(J) is different from the as-observed 
action distribution. The fraction of its orbital period that a 
given star spends in the Solar neighbonrhood depends on its 
orbit, so that the probability of finding the star in the Solar 
neighbourhood depends on its action. In section [T] we take 
into account this probability to reconstruct the observable 
distribution function /q(J). 


3.2 Observable vertical velocity dispersion 

In order to have a physical insight into T©, let us consider 
the velocity space that corresponds to Tq by assuming a po¬ 
tential model of the Milky Way. Figure [4] shows the allowed 
range of [vr, v^) for disc stars located at {R, z) = (6, 0) kpc 
with zero vertical velocity to reach the Solar neighbourhood. 
Since only a portion of the disc stars at 7? = 6 kpc can reach 
the Solar neighbourhood, what we can learn from /© about 
the thick disc at i? = 6 kpc is the velocity distribution of 
stars within the ‘observable’ region in Figure 21 

With the information on /q, we can estimate the ob¬ 


servable vertical velocity dispersion defined by 


o'z,0(-R,O) 


fd^v fQ{J) 


( 2 ) 


at any given radius R on the disc plane. This dispersion 
corresponds to the vertical velocity dispersion of those stars 
located at R whose velocities are within the observable re¬ 
gion. Thus Uz.Q is conceptually different from the (ordinary) 
vertical velocity dispersion given by 


al{R,0) 


f d^v vlf{J) 
jd^v fiJ) 


( 3 ) 


which is the vertical velocity dispersion of all the stars lo¬ 
cated at R. 


4 METHOD 

Here we describe how to reconstruct /©(T) of the thick 
disc from a nearby sample of thic k -disc stars, based on 
the idea of Sommer-La rsen fc ZherJ (Il990l b We note that 
iMcMillan fc BinnevI 1 2 OI 2 I 1 investigated a rigorous treat¬ 
ment of an as-observed distribution function in which survey 
selection function is taken into account (see Appendix). We 
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Figure 4. The allowed region in the (n/?, n^)-space for disc stars 
at {R, z) = (6, 0) kpc with = 0 km s~^ to reach the Solar neigh¬ 
bourhood, calculated for the potential model in[3 The boundary 
curve between the observable region and unobservable region cor¬ 
responds to -Rapo = Ro = 8-3 kpc (Rapo is the apocentric distance 
of the orbit), and thus the circular orbit at R = 6 kpc (red dot) 
lies within the unobservable region. When \vz\ > 0 kms“^, the 
boundary curve slightly shifts downwards in this diagram. How¬ 
ever, this difference is not significant for a typical value of \vz \ of 
the thick-disc stars. In this plot, we assume a realistic potential 
model of the Milky Way described in section [5] 


By following ISommer-Larsen &; ZlienI (|l990l ). we ap¬ 
proximate the observable distribution function /q to be the 
weighted sum of single orbit distribution functions given by 

N 

feiJ) ( 6 ) 

1=1 

Here, Si is a selection index that determines whether the 
ith star should be included (si = 1) or not (si = 0), 
which shall be introduced to stabilise our analyses. The 
factor Wi{> 0) i s the orbit weighting fa c tor to be deter¬ 
mined. Although ISommer-Larsen fc ZhenI (ll99Clfl evaluated 
Wi through Bayesian statistics, here we numerically evaluate 
p{Ji\Xi) and simply set Wi = 1/p[Ji-,Xi). This simple choice 
of Wi corresponds to a limiting ca se where all of our sam¬ 
ple s tars have identical position ilSommer-Larsen fc ZhenI 
I1990Ii FI Since most of our sample stars are located within 
100 pc away from the Sun, our simple choice of Wi is well 
justified. 

4.2 Reconstruction of vertical velocity dispersion 
of the thick disc 

Once /q is reconstructed, we can estimate the observable 
vertical velocity dispersion Uz,© at a given point on the 
disc plane, by using equation ([2]). It can be shown that 
^■^.©(RoiO) = (Jz(Ro,0) if the assumed potential is correct 
and a large enough number of local sample stars are avail¬ 
able. Thus, we can perform a sanity check of the adopted 
potential by seeing if these dispersions reasonably agree with 
each other. 


do not follow their formulation here since we do not know 
the exact survey selection function for our sample stars (but 
we do know that our sample stars are not kinematically se¬ 
lected). Although we do not take into account the survey 
selection function, our approach has an advantage that we 
do not have to assume the functional form of the distribution 
function. 


4.1 Reconstruction of observable distribution 
function /© of the thick disc 


In the following, we regard disc stars as test particles moving 
in the Milky Way potential. Also, we assume that we have 
a kinematically-unbiased sample of N thick-disc stars which 
happen to be passing near the Sun. 

Under an axisymmetric pote ntial 4>(R, 2:), each star has 
three actions J = {Jr, J^, Jz) (iBinnev fc Tremaine! |2004 
|Binnevl [201^). Let us denote the actions of the ith star as 
Ji. Then, the time-averaged distribution function of a single 
orbit characterised by Ji is expressed as 


fr^^{x,v) 


-^S{J{x,v) - Ji), 
(2tv) 


and the corresponding spatial density is given by 


p{Ji-,x) = J ffv 


1 

dv 

(27r)3 

dJ 


( 4 ) 

( 5 ) 


The value of p{Ji;x) can be interpreted as the probability 
density of finding a s tar with action Jj at position x, since 
/ (fxdPv = 1 dchiba fc BeerdlioOlli . 


5 ANALYSIS AND RESULTS 


5.1 Analysis: Derivation of weighting factors 


We adopt a realistic potential model of the Milky Way, con¬ 
sisting of halo, bulge, and disc components. 

In our model, the h alo potential is given by the Navarro- 
Frenk-White potential (iNavarro et al.l[l996lj of the form 


2 In (1 -I- r dn) 

4«halo = -Vii - 

r/dR 


( 7 ) 


where r is the Galactocentric distance, uh = 379.79 km s ^ 
and dn = 9.1 kpc. The bulge potential is given by 


4*bulge — 


GMb 

cb +r 


(8) 


dHernauisd Il99fll l. We adopt Mb = 2.05 x lO^^M© and 
Cb = 2 kpc. We have confirmed that adopting a slightly 
flattened bulge potential hardly affects our results. For the 
disc potential, we assum e the Miyamoto-Nagai potential 
llMivamoto fc Na^ll975l ~l given by 


4-h 


GMo 


\jdJJ + (Rd + -I- Zp) 


( 9 ) 


^ This property can be intuitively understood in the following 
manner. The probability p of finding a given star in the Solar 
neighbourhood is equal to the fraction of its orbital period to 
stay in the Solar neighbourhood. Thus, if such a star is observed 
in the Solar neighbourhood, it represents 1/p stars in the same 
orbit (including those with different orbital phases). 
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where Md = 5.77 x 10 ^'^Mq, Rd = 3.2 kpc, and zu = 
0.23 kpc. The values of parameters are chosen by using the 
publicly availabl43 code of galpy (IBov ^ l2015ll . while fixing 
Ro = 8.3 kpc, ulsr = 240 km s“^, and cb = 2 kpc. 

After setting up the potential models, we derived the 
weighting factor for the ith star, Wi, in the following manner. 
We first integrated the orbit of the ith star forward in time 
for a long enough period. We then calculated the fraction 
Sf of the integration time which the star spent in a volume 
defined by Ri — SR < R < Ri + SR, Zi — Sz < z < Zi + Sz, 
and 0 ^ (/) < 27r. Here, we denote the observed position 
of the ith star on the meridional plane as {Ri,Zi) and we 
adopted {SR,Sz) = (0.05,0.1) kpc. By taking into account 
the relation Sf = p{Ii', Xi) x inRiSRSz, we derived p(/i; Xi) 
and set Wi = Xjp(Ii\ Xi). We set Si = 0 for those stars that 
satisfy Wi > ({uij})+3a'{uj3}, where ({wj}) and are the 
mean and the standard deviation of {wj}, respectively. This 
final procedure was included to prevent our results being 
affected by small number of stars with large Wi. 


5.2 Reconstructed vertical velocity dispersion 


Figure [5] shows the reconstructed profile of (Jz,q- Here we 
use an approximation of 


Ei SiWi dz dR' 2nR' J d^v vlff 


0-^,0 ( i ?, 0 ) ~ 


Ei SiWi dz dR' 2nR' f d^v f‘ 


single 


( 10 ) 


We adopt {AR,Az) = (0.25,0.1) kpc and calculate 

az^ 0 {R,O) at 4.5 ^ 7?/kpc 9.5 in steps of 0.5 kpc. 
The radial range of R is set so that more than 42 stars 
(i.e., more than N/3 stars) have non-zero contribution to 
the integrals in equation IIOJ, i,e., have positive values of 
2nR'Jd^vfr^^^). 

In Figure [5l the red dots show the reconstructed values 
of (72,0(7?, 0), and the shaded region represents the uncer¬ 
tainty which is derived from a bootstrap resampling of the 
contributing stars. The velocity dispersion declines nearly 
exponentially as a function of R, and the best fit exponen¬ 
tial profile of (72,0 is shown with a black straight line. We 
find that the best-fit scale length defined by 

( 11 ) 

is Rs = 8.3 ± l.l(rand.) ± 1.6(sys.) kpc. The random er¬ 
ror is the uncertainty in the fit and the systematic error is 
estimated by performing mock catalog analyses. 

We interpolate the values of ctz.q ai R = 8.0 kpc and 
R = 8.5 kpc to obtain CT2,0(7 ?o,O) = 39.4 ± 2.5 kms“^. 
This value is in agreement with the local value of — 
41.6 km s“^ in our sample, indicating that the adopted po¬ 
tential is reasonable at least in the Solar neighbourhood. 



R [kpc] 

Figure 5. The reconstructed observable vertical velocity disper¬ 
sion (72,0 of the thick disc as a function of R. 


nearly exponential profile of (Jz.q is a generic feature inde¬ 
pendent of the assumed parameters. Rs mildly depends on 
these parameters, but it only increases from 7.5 kpc to 9.4 
kpc when we vary Ro from 2 kpc to 4 kpc; and from 8.1 kpc 
to 9.1 kpc when we vary zo from 0.15 kpc to 0.5 kpc. 


5.3 Dependence of the vertical kinematics on vr 
and 

As discussed in section [331 o'2,0(7?, 0) represents the verti¬ 
cal velocity dispersion of stars located at {R, 0) whose ve¬ 
locity is within the observable region in the velocity space 
(Figure Ul. In order to see the variation of the vertical kine¬ 
matics within the observable region, we create the follow¬ 
ing groups of stars from our sample: (1) those stars with 
5.5 < i?peri/kpc < 6.5 {N = 26 stars); and (2) those stars 
with 5.5 < (l/2)(i?peri + 7iapo)/kpc < 6.5 (N = 30 stars). 
Here, we define Rperi and i?apo as the minimum and maxi¬ 
mum value of R for each star’s orbit. These samples happen 
to be non-overlapping, and the observed vertical velocity dis¬ 
persion for groups 1 and 2 are respectively 38.9±5.4 km s“^ 
and 42.6 ± 5.5 km s“^. All of these stars spend some time 
at {R,z) = (6.5,0) kpc, so in Figure |6l( a), we show the dis¬ 
tribution of I Ur I and when they pass through this loca¬ 
tion. The stars in groups 1 and 2 are distributed in distinct 
areas in this velocity space, reflecting the different orbital 
eccentricities. In Figure |5Kb), we show the vertical velocity 
dispersion CT2,0(7?, 0) of these groups, calculated by using 
{uii} derived in section [5] We see that the radial profiles of 
(72,0(7?, 0) for these groups are quite similar to each other 
at 6 ^ 7?/kpc ^ 8.5. This result suggests that the vertical 
kinematics of the stars in the observable region is more or 
less homogeneous at least at 6 ^ 7?/kpc ^ 8.5. 


5.2.1 Dependency on the adopted potential parameters 

We vary the parameters (7?d,zd) within the range of 2 < 
7?D/kpc < 4 and 0.15 < 2 d /kpc < 0.5 and find that the 

® Available at http://github.com/jobovy/galpy 


6 DISCUSSION 

6.1 Dependency of the vertical kinematics on 
radial and azimuthal velocities 

As seen in Figures [21(a), (21^), and [3] the vertical velocity 
dispersion of the thick disc does not have a strong depen¬ 
dency on the in-plane velocities [vr and v^,) in the Solar 
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Figure 7. The kinematically-estimated thickness of the thick 
disc as a function of R by assuming az{R, 0) = az,Q{R, 0). Also 
shown is the kinematically-estimated thickness of the local thin 
disc which is calculated by assuming = 20 it 2.5 km s“^. 



Figure 6. (a) The distribution of stars in group 1 (5.5 < 

^peri/kpc < 6.5) and group 2 (5.5 < (l/2)(ijperi + Rapo)/kpc < 
6.5) in the (Ui^Jui}!) plane when they pass through {R, z) = 
(6.5, 0) kpc. (b) The reconstructed profile of cr^ 0 for these groups. 
As a reference, we show the exponential profile shown in Figure 
[ 5 ] (black solid line). 


neighbourhood. In order to have some physical insight into 
this finding, let us suppose, as a working hypothesis, that 
this apparent homogeneity of vertical velocity distribution 
across the (wk, u,^)-space is the case not only near the Sun 
but also outside the Solar vicinity. Then it is implied that 
the distribution of the vertical action Jz is not strongly 
dependent on the radial action [Jr) or azimuthal action 
{J 4 >)- It follows that there is a typical value of vertical ac¬ 
tion Jz,typical that characterises the vertical motions of the 
thick-disc stars. On the other hand, the maximum height 
Zmax above the disc plane that a given thick-disc star can 
reach at Galactocentric radius R is determined mainly by its 
vertical action Jz ■ In realistic potential models of the Milky 
Way, Zmax — Zmax{Jz,R) is an increasing function of R, 
since the gravitational pull towards the disc plane is weaker 
at larger R. Thus, if the distribution of Jz of the thick-disc 
stars is nearly independent of {Jr, JV), then the thick disc 


needs to flare and the scale height of the thick disc roughly 
follows the radial behaviour of Zmax{Jz,typical, R)- 

Then, under what circumstances does the distri¬ 
bution of Jz become (approximately) independent of 
{Jr, J<i>V Here we consider such an example. Let us 
suppose th at the thick disc was fo r med through radia l 
migratio n (ISellwood &: BinnevI 20021: Roskar et al.l l2008l : 


ISchonrich fc Binnev 2009a Loebman et al.l 201ll)~ Some 
numerical simulations suggest that the radial migration con¬ 
serves the stella r vertical action Jz as the star changes its 
in-plane motion JSolwav et akllioi^ ). Thus, radial migration 
redistributes the stars in {Jr, J,^)-space while keeping Jz un¬ 
changed. This redistribution in the action space is expected 
to erase the initial dep endency of Jz-distribution on {Jr, J,/,) 
dMinchev et al.|[2012l ). unless the initial dependency is too 
strong to erase. If we take into acco unt that the thick-dis c 
stars are generally very old (> 8 Gvr: lHavwood et al.ll201^ ). 
it may be likely that the radial migration ends up with a 
homogeneous distribution of Jz as a function of Jr or J^. 
It is of interest to note that some numerical simulations 
predict that the thick disc flares if it was formed though 
radial miEration (ILoebman et al.ll201ll:lMinchev et al.ll2012l; 
iRoskaretabim^ 

If (Tz is really independent of {vR,vf) anywhere on the 
disc plane, we have Gz{R,0) — (Jz,q{R,0)- In this case we 
can estimate the thickness H{R) of the thick disc at Galac¬ 
tocentric radius R by using our reconstructed profile of ctz.q. 
As an example, we define the thickness H by 


^{R,H{R)) = ^{R,Q) + -ol{R,Q), 


( 12 ) 


where $(i?, z) is the Galactic potential. By plugging Gz,q 
into CTz in equation inj, we obtain H{R) shown in Figured 
Although the H{R) profile mildly depends on the adopted 
potential, we have confirmed that H increases as a function 
of R independent of the potential parameters as long as they 
are within the range explored in section [5.2.1l 


6.2 How can we keep the scale height constant? 

As we have discussed in the previous subsection, the radial 
dependence of the thickness of the thick disc is related to 
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the dependency of the vertical kinematics on the in-plane 
motions. 

By taking into account some observational evidence 
of the thick-disc compone nt in external galaxie s with 
a constant scale height (Ivan der Kr uit fc Searld Il98ll : 
Ijensen &: Thuarilll982l : IComeron et al.ll2011^ . let us consider 
what we would need to keep the scale height H constant. 
Here we concentrate our discussion on the situation at a cer¬ 
tain radius R = Ri in the inner disc {R\ < Rq). In order to 
satisfy H{Ri) = H{Ro), we require az{Ri,0) > ctj, 0(J?i,O) 
so that the vertical velocity distribution is inhomogeneous in 
the {vr, v^) space. Since crz,Q represents the vertical velocity 
distribution of observable stars (those stars whose velocity 
is inside the observable region) and Oz represents the ver¬ 
tical velocity distribution of all the stars, this requirement 
implies the following situation: In the inner disc, those thick- 
disc stars with large J 4 , (large guiding centre radii) need to 
have smaller vertical action than those thick-disc stars with 
small and This situation may be attained if there are 
enough low-eccentricity, shell-like orbits with large vertical 
motions in the inner disc. 


7 CONCLUSIONS 

We have investigated an orbit-based method to reconstruct 
the observable distribution function /q , which describes the 
distribution of action within the portion of action space cov¬ 
ered by a local disc sample. By applying our method to 127 
chemically-selected thick-disc stars in the Solar neighbour¬ 
hood, we found that the vertical velocity dispersion (Jz,q 
that corresponds to /q decreases nearly exponentially as a 
function of R with a scale length of Rs = 8.3 ± l.l(rand.) ± 
1.6(sys.) kpc. 

We also found that the vertical velocity dispersion Uz 
of the local thick-disc stars shows only weak dependency on 
the in-plane motions [vr, Vrj,). This apparent homogeneity of 
(Tz in {vr,v^)- spa.ee may be naturally explained if the dis¬ 
tribution of the vertical action Jz is only weakly dependent 
on the radial action Jr and the azimuthal angular momen¬ 
tum Jcj,. One possible mechanism to produce such action 
distribution is radial migration, where the thick-disc stars 
are re-distribu ted in the (Jr, J<<,)-sp ace while keeping their 
Jz unchanged llMinchev et al.ir2012f l. 
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APPENDIX A: SELECTION FUNCTION OF 
THE SURVEY 

In our analyses, w e foll owed the approach of 
ISommer-Larsen & ZhenI (ll99Cll l to relate the observed 
kinematic data to the underlying distribution function. 
We note here that our approach owes much to the fact 
that our sample stars are located at the immediate Solar 
neighbourhood (d < 100 pc). 

Let us consider a sample of stars taken from an imag¬ 
inary magnitude-limited survey in which all the stars with 
apparent magnitude brighter than a certain value are ob¬ 
served. In this case, those sample stars at large distances 
are high-luminosity stars, while the sample covers a wider 
luminosity range for nearby stars. Thus, such a sample puts 
more weight on nearby stars. If the kinematical properties 
are highly inhomogeneous within the surveyed region, this 
sample is kinematically biased. In general, if the surveyed 
region is large, the survey selection function has to be taken 
into account in order to estimate the kinematic properties 
of the system (McMillan & Binney 2012). This means that 
our methodology is not applicable to a sample that covers a 
large volume unless the completeness is taken into account, 
since our method implicitly requires that the sample stars 
are fair representatives of the stellar population. However, 
since our sample is distributed in a small region (d < 100 
pc), it is safe to assume that the kinematical properties of 
the thick-disc stars are homogeneous inside this region (the 
thick disc has a scale height of ~ 1 kpc and a scale length 
of 2-3 kpc). Therefore, the effect from the (poorly-known) 
survey selection function is relatively small in our sample. 
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